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Abstract. An axion with a mass greater than 1 eV should be detectable through its 
decay into two photons. In this paper I discuss the astrophysical and cosmological limits 
which define a small window of allowed axion mass above 3eV. A firm upper bound to 
the axion’s mass of m a < 8 eV is derived by considering the effect of decaying axions 
upon the diffuse extragalactic background radiation and the brightness of the night sky 
due to axions in the halo of our galaxy. The intergalactic light of clusters of galaxies is 
shown to be an ideal place to search for an emission line arising from the radiative decay 
of axions. An unsuccessful search for this emission line in three clusters of galaxies is 
then detailed. Limits to the presence of any intracluster line emission are derived with the 
result that axions with masses between 3eV and 8eV are excluded by the data, effectively 
closing this window of axion mass, unless a severe cancellation of axionic decay amplitudes 
occurs. The intracluster flux limits are then used to constrain the amplitude of any such 
model dependence. 
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I. Introduction 

Astrophysical and cosmological considerations have lead to numerous constraints on 
new particle physics phenomena. Nowhere has the interplay of these disciplines been more 
apparent than in the search for the axion. 1,2 The axion is a direct consequence of a modest 
and desirable extension of the standard model of particle physics, Peccei-Quinn (PQ) 
symmetry. 3 PQ sy;~metry was originally proposed as a solution to the string CP problem 
of Quantum Chromodynamics (QCD). PQ symmetry also arises naturally in theories of 
supersymmetry and superstrings. 

As a result of non-trivial topology of the QCD vacuum, the strong interactions violate 
CP invariance with a strength proportional to a parameter ©. This CP violation manifests 
itself in a n um ber of ways, most notably in the electric dipole moment of the neutron, 
d n . QCD predicts a value of roughly d n ~ 5 x lO -16 0e cm while the experimentally 
determined value is d n < 10 -25 e cm. This leads to the constraint 0 < 10 -9 — 10" 10 . A 
value of 0 this small violates the naive expectation that a dimensionless parameter of a 
theory should have a value 0(1) and is unnatural in the sense defined by ‘tHooft. 4 ’ 6 The 
strong CP problem is, in essence, the question: why is 0 so small? 

The PQ solution solves this problem by introducing a new global U( 1) symmetry into 
the theory. This 17(1 )pq is then spontaneously broken at an energy scale /pq, where I 
follow the notation of reference 1. The pseudo-Nambu-Goldstone boson associated with 
this sy mm etry breaking is the axion. 6 U{ 1)pq is not an exact symmetry; it is broken by 
QCD ano malies at the quantum level. Because of this explicit breaking the axion acquires 
a small mass which is inversely proportional to the scale qf spontaneous symmetry breaking 

sfz f n m„ 0.62 eV • 10 7 GeV 

= — = — 7^n — • (1) 

where N is the color anomaly of the PQ symmetry, 2 = m u /rrii ~ 0.56, /*• = 93 MeV , 
and m„ = 135 MeV . (Throughout I refer to the axion’s mass in terms of its rest energy, 
m a c 2 .) In addition, the axion’s couplings to matter axe inversely proportional to /pq (or 
equivalently, proportional to m a ). At the time of the QCD confinement phase transition 
the axion potential develops a minimum which corresponds to 0 = 0. As the axion relaxes 
to the mi nimum of its potential it forces 0 = 0, thus solving the strong CP problem. 

There are two generic types of axions. DFSZ axions couple to both quarks and leptons 
at tree level. 7 The hadronic, or KSVZ, axion has no tree level couplings to leptons but 
does couple to them through loop diagrams. 8 Throughout the remainder of this paper I 
will focus on the hadronic axion, as DFSZ axions more massive than 1 eV have been ruled 
out by stellar evolution arguments. 1,2 Since DFSZ axions couple directly to electrons, they 
can be profusely emitted from the cores of stars. This provides a very efficient means of 
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transporting energy from the core to interstellar space. For masses greater than 1 eV a 
DFSZ axion would so strongly effect the core of the Sun that its existence could not have 
gone unnoticed. This lower bound to m a (DFSZ) is further strengthened by considering 
the effect of axions upon red giant stars. These mass bounds on DFSZ axions escape much 
of the model dependence which is present in similar bounds to the hadronic axion (as 
discussed in the next se^ion) and are much stronger since the direct aee coupling is much 
stronger than the equivalent coupling for the hadronic axion. 

Axions can decay into two photons. The 077 coupling arises due to two different decay 
mechanisms: through axion-pion mixing and via the electromagnetic (EM) anomaly of PQ 
symmetry 
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— J& ww " u * )dS - B 

where a ~ 1/137, F* u is the EM field strength tensor and jF m „ is its dual, 1.95 = 2(4 + 
z)/3(l + z), and E is the value of the EM anomaly of PQ symmetry. The axion decay 
width to two photons is 


r = t -1 = 

CL ' n 
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2 s tt 3 ( fpq/Ny 


( E/N - 1.95) 2 ml 


which leads to an axion lifetime of 


r a (a — ► 27) ~ 6.8 x 10 24 ( 2 (m a / eV) 5 sec (2) 

where I have introduced the notation ( = (E/N — 1.95)/0.72. In the simplest models which 
incorporate axions, E/N = 8/3 and C = 1. While E/N = 8/3 is the value that arises when 
the axion is incorporated into the simplest unified models it is not required; models with, 
for example, E/N = 2 are easily constructed. 9 For models of this type, C = °- 07 the 
two photon decay of the axion is highly suppressed due to the near cancellation of the 
amplitudes from the EM anomaly and axion-pion mixing. The parameter ( is of great 
importance when considering astrophysical limits to the mass of the axion. 

While the properties of the axion are well defined in terms of /pq (or m a ), this quan- 
tity has no preferred value. / PQ might plausibly lie anywhere in the range 10 2 GeV to 
10 19 GeV, corresponding to axion masses of 100 keV <m,< 10 -12 eV. Astrophysical and 
cosmological arguments have been extremely useful in narrowing this allowed range. Due 
to many clever arguments there are but two “windows” of mass left open to the axion: 
10 -6 eV < m a < 10~ 3 eV and, for the hadronic axion only, 3eV < m a < 8eV. The latter 
window is the subject of this paper. In the next section I will review the arguments that 
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define the boundaries of the multi-eV window; more complete reviews appear in references 
1 and 2. 

Axions, if they exist in either of the windows of allowed mass, will have cosmologi- 
cally interesting and potentially detectable relic abundances. The dominant production 
mechanisms in the two windows differ significantly. In the lower mass window, axions 
are produced through two different coherent mechanisms: misalignment production and 
axionic string decay. Both of these processes are highly non-thermal. Axion production 
via the misalignment mechanism produces axions with an abundance 10,11 

- 0.7 


ft a (mis)/i 5 0 ~ 0.52 x 10 


±o.4 ( A QCD \ 
\ 200 MeV ) 


<*> 


where 9\ is the initial angular displacement, between 0 and 7r, of the axion field from its 
mi nimum, 0 = 0, A qcd — 150-400 MeV is the QCD scale parameter, and f{6\) is an 
increasing function of Q\ (/( 0) = 1). If an inflationary epoch occurred after the PQ phase 
transition, 6i will have a unique value throughout the observable Universe and cc Q\. If 
there was no inflation, or the PQ transition occurred after the inflationary epoch, Q\ will 
have a different value in each causally distinct region of the Universe and the correct value 
of 9i to use in eq. (3a) is the rms, 0\ = irjs/Z. Inserting this into eq. (3a), one finds 


f2 a (mis)^ 0 ~ 3.4 x 10 ±0 ' 4 


( Aqcd \ °' 7 / rUg \ 
V 200 MeV ) \10 -s eV/ 


- 1.18 


(36) 


If axions are dominately produced via axionic string decay (which cannot happen if the 
Universe inflates after the PQ transition) their abundance is 


fl a (string)/i 5 0 ~ 



m a \ 
10“ 6 eV / 


- 1.18 


or 


/ m a \-ii8 
\10 -3 eV / 


(3c) 


There is still considerable controversy as to which value is correct. (For a more thorough 
discussion of these axion production mechanisms see references 1, 2, and 10.) By examining 
Eqs. (3a-c) it is apparent that if the axion exists with a mass in the 10 -6 eV-10~ 3 eV 
window, it likely will contribute a fractional density of ~ 0(1)- Thus an axion in this 
window will make up a sizable fraction of the dark matter in the Universe. Axions with 
masses above 10 -2 eV are produced thermally. Since the coupling of the axion to normal 
matter is proportional to m a , axions eventually couple strongly enough to thermalize. 
This happens for axions with masses greater than 10 -2 eV. 10,12 Thus, axions in the multi- 
eV window are thermally produced and have an abundance similar to that of microwave 
background photons and relic neutrinos. 

Once these axions are produced, where are they likely to be found? Since the current 
temperature of the microwave background radiation is T 7 o — 2.4 x 10 -4 eV, axions are 
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non-relativistic and have been since before decoupling. Therefore axions should, in accord 
with the equivalence principle, fall with baryons and any other particles into the various 
potential wells which develop in the Universe. The most likely place to find relic axions 
with masses of a few eV is in clusters of galaxies and the halos of galaxies (as the axions 
have no way of dissipating their energy to condense further). As I will show later, phase 
space considerations make clusters of galaxies the most likely hunting ground for th^e 
axions. 

If these multi-eV axions reside in clusters of galaxies, as seems likely, their decay 
into two photons provides a powerful way of searching for them. 12 ’ 13 Decaying axions 
will produce an emission line at a wavelength A„(z) ~ 24800 A(eV/m a )(l + z) where 
z is the cluster’s red shift. This line will be Doppler broadened by the velocities ax- 
ions have in the cluster to a width of AA a ~ 100 A and will have an intensity of I a ~ 

I. 5 x 10 - 17 (m a /3eV ) 7 £ 2 erg cm -2 arcsec -2 A -1 s _1 . In any attempt to detect this feature 
the signal must compete with the brightness of the night sky (NS). The night sky is charac- 
terized by a continuum intensity level of 2ns ~ 10 -17 erg cm -2 arcsec 2 A 1 s 1 and many 
strong atmospheric emission lines 14 (see Figs. 2-6). In this paper I give a detailed account 
of an unsuccessful telescopic search for an axionic decay line in three clusters of galaxies. 
The search was conducted at Kitt Peak National Observatory (KPNO) and no axion decay 
1in<» was detected, effectively closing the 3 eV to 8 eV window of axion mass. 

I close the introductory comments by noting that while the search was motivated by 
the axion, our results are quite general. They can place meaningful constraints upon any 
relic particle species which decays to one or more photons and which clusters with galaxies, 
e.g., massive neutrinos. 15 Even more generally, the results place limits on the presence of 
any diffuse inter-galactic line emission in clusters of galaxies, no matter what the source. 

II. The Window 

The 3 eV to 8 eV axion window is defined by a number of different arguments. The 
lower mass boundary is determined by the axion’s effect upon SN 1987A and its observed 
neutrino burst. 16 The upper boundary is determined by two very different considerations. 
The first of these is the effect of axion emission upon the life cycle of red-giant (RG) and 
helium burning (HB) stars. 17 The second is the effect of decaying axions upon the diffuse 

extragalactic background radiation (DEBRA). 12,19 

The effects of axions upon SN 1987 A have been considered by numerous authors 
(see references 1,2, and references therein) with the result that axions in the mass range 
10 -3 eV < m a < 3 eV are not compatible with the observations and therefore excluded. 
Here I will sketch the argument which leads to the upper bound of 3 eV (for details see 
reference 16). As the axion’s mass is increased, it couples more and more strongly to “nor- 
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mal” matter. At masses greater than about 10 -3 eV the axion couples strongly enough so 
that it is profusely emitted by the core of a hot young neutron star during a supernova 
explosion, but no so strongly so that the matter seriously impedes the axion energy flux 
out of the star. This flow of energy out of the core severely depletes the resevoir of heat 
which powers the later phase of neutrino emission, severely shortening the neutrino burst. 
As the axion mass is increased beyc.id 10 -3 eV the axion couples even more strongly and 
at a mass of about 10 -2 eV the axion’s mean free path becomes less than the radius of the 
neutron star; axions axe trapped, much like the neutrinos. This axion trapping reduces 
the flux of axion energy leaving the neutron star. As the axion mass is further increased 
the trapping becomes so strong that, like photons, the axions carry too little energy to 
seriously effect the neutrino burst. 

The greatest observable effect of axions upon the supernova neutrino burst is in the 
burst duration. By incorporating axion emission into numerical models of proto-neutron 
star formation, the authors of ref. 16 determined that the SN 1987A neutrino burst 
observed in the KII and IMB detectors would be shortened by a factor 2 or more for 
axions with masses less than 3 eV (and more than 10 -3 eV). This constitutes the lower 
bound of the multi-eV window. There axe a number of uncertainties in this limit. For 
example, the exclusion criterion of requiring the neutrino burst to decrease by a factor of 
two, while reasonable, is rather arbitrary. Also, while the relevant axion-nucleon couplings 
are relatively model independent they are undetermined to factors of 0(1). Taking these 
and other concerns into account the authors of ref. 16 estimate an uncertainty of a factor 
2 in their result, meaning the true lower bound to the window lies in the range 1.5 eV to 
6 eV, with 3 eV being the probable value. 

Stellar evolution arguments provides one means of defining the upper bound to the 
multi-eV axion window. Raffelt and Dearborn 17 have considered the effects of hadronic 
axions upon various phases of a star’s evolution (see also refs. 18 and 19). In their paper 
they derive two limits to the mass of the axion which are relevant to the present search. 
The first is m a < 2£ -1 eV from considering the effects of axion emission upon the peak 
brightness of RG stars at the tip of the giant branch. The second limit is even more 
restrictive, m a < 0.7£ _1 eV. This constraint derives from the fact that axions remove 
energy from the core of HB stars at a sufficient rate to shorten their calculated ages by a 
factor of 2. This would cause a discrepancy with the observed number of HB stars in the 
open star cluster M67. Both of these constraints seem to rule out the claimed existence of 
the multi-eV window (when combined with the results from the SN 1987A analysis) but 
there are at least two reasons to consider the window further. 

Both of the above constraints depend upon the Primakoff effect for axion emission 
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from the stars. This effect makes use of the same anomaly diagram which leads to the two 
photon decay of the axion. Hence the same model dependence, as embodied by (, which 
enters the axion lifetime plays a role here. In fact letting E/N = 2 (C = 0.07), the above 
limits become m a < 30 eV and m a < 10 eV respectively; well in excess of the lower bound 
from SN 1987A. We will find that we can do much better than this from considering relic 
decays. The other reason to consider the window further is simply that both of the above 
arguments are based upon the statistics of small numbers. Unlike the bound from SN 
1987A, where the 19 observed neutrino events over about a dozen seconds trace the entire 
cooling history of the proto-neutron stax, RG and HB stars live for many millions of years. 
Thus, in order to infer the properties of RG and HB stars one must argue based upon the 
observed number of such objects as compared to other types of stars. In both of the above 
cases the data is a small sample which may, or may not, be statistically significant (in M67 
there are only 5 HB stars). So, while the physics upon which the above two bounds are 
based seems quite sound, the observations are not yet at a level to confirm them. Given 
the above possible loopholes and uncertainties in the stellar evolution arguments it makes 
sense to consider what other upper bounds one may derive to the mass of the axion. On 
possible method is to consider the effects of decaying axions upon the DEBRA. 

Measurements of the DEBRA flux have, in several instances, been used to constrain the 
properties of neutrinos. 20 ’ 22 In much the same way these flux measurements can be used 
to constrain the properties of the axion. 19 If one makes the most conservative assumption, 
that axions are distributed uniformly throughout space, axions produce a diffuse glow in 
the night sky with an intensity of 12,22 


L = 


d? E 


n a m a c 3 (A a /A) 3 


n a m a c 


7/2 

^ ) ( 4 «) 


dAdQdXdt 47rifor a A a [l — ^tot + ^TOT(A/A a )] x / 2 4irHoT a X 0 
valid for A > A a , and where n a = present axion number density, H 0 = the Hubble constant, 
A = observed wavelength, A a = 24800 A( eV/m a ), and Otot = P/Pcrit is the universal 
density parameter. The second step in eq. (3a) follows from setting ft T OT = 1, the 
theoretically preferred value. It is interesting to note that if an axion with a mass of 
6 eV or so was uniformly distributed throughout the cosmos, the spectral shape of the 
resulting glow might present a powerful method of determining Qtot- Unfortunately, as 
I will show, there is no such feature. From now on I will assume fi T OT = 1 and use the 
second expression in eq. (3a). 

If I now insert values of the quantities into eq. (3a) I find that an axion will produce 
a glow in the night sky of 12 

7/2 


I a ( A > A a ) = 1.8 x 10 


-23 ( 


/m a > 

I 7 C ' 1 

(}L\ (b±\ 

V^v> 

’ /*50 



erg cm 


' 2 arcsec 2 A 1 s 1 (46) 
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where g»f is the number of relativistic degrees of freedom when the axion freezes out (see 
the next section). Note that any limit to the axion mass derived from this will depend on 
C -2 / 7 ; this is a much weaker dependence than the stellar evolution limits. It is now a simple 
matter to compare the flux predicted in eq.(3) to that measured in DEBRA experiments. 

There axe numerous experiments which have searched for and/or found the intensity 
of the DEBRA. For a recent review of the available data see ref. 20. I will focus on 
results which can be derived for axions with masses between 3 eV and 25eV, the size of 
the window if we neglect the stellar evolution arguments. There are a number of reasons 
for restricting ourselves to this region. First, SN 1987A precludes the existence of axions 
with m a < 3 eV. Second, as I will show, no useful constraints on m a may be derived from 
the DEBRA for m a < 3eV. Third, for m a > 25 eV, r 0 is shorter than the age of the 
Universe and eq. (4) is no longer valid. Lastly, for axion masses greater than about 25 eV 
the wavelength of the axion decay radiation is shorter than the Lyman limit, and a truly 
extra-galactic flux cannot be measured due to the opacity of the interstellar medium. In 
this range I may now compare the intensity predicted by eq. (4) to that actually measured. 
In Fig. 1 I present results of this exercise. I consider the intensity, I a , for a number of 
axion masses and require a value of £ small enough so that I a ( X = A a ) = /debra- It is 
clear that throughout the region of interest £ must be very small to be consistent with the 
observations. In fact, £ = 0.07 ( E/N = 2) occurs at an axion mass of 7.5 eV. This value 
of m 0 constitutes the upper bound to the multi-eV window of axion mass. Fig. 1 also 
clearly shows that as m a approaches 3 eV, £ goes to 1. The lowest value of m a for which 
the available data gives any information is m a = 3.8 eV. For this mass £ < 0.43, so that 
in the simplest axion models consideration of the DEBRA flux requires m a < 3.8 eV, this 
is still in excess of the SN 1987A bound. 

The expected signal from axions residing in the halo of our galaxy provides a compli- 
mentary mass bound to that deriving from the DEBRA flux. This is because these two 
limits span a wide range of conditions — from very clustered to totally unclustered. The 
galaxy analysis will thus strengthen the upper bound to the window of m a < 8eV. 

The density of the galactic halo is usually parameterized as 


P(r) = P© 


R% + a 2 

r 2 + a 2 


(5) 


where p© is the local halo density, i?© = 8.5 kpc is the distance from the Sun to the 
galactic center, and a is the halo core radius. Based upon galactic modeling, p© ~ 5 x 
10 -25 gcm -3 ~ 0.0lM©pc -3 and the halo core radius lies between 2 kpc and 8 kpc. 23 Note 
that, for r a, the density run of eq. (5) p(r) oc r -2 , just as one would expect based 
upon the observed flat rotation curve of our galaxy. The surface density observed in the 
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model of eq. (5) is E(M) = R®p®J(x,9) where x = g/Rq, 9 is the angle between the 
galactic center and the line of sight, and 


J(x,0) = 


1 + x 2 


y/ x 2 + sin 2 0 \ ^ 


7T _ i 

— “(~ tan 


cos 9 


a/x 2 + sin 2 0 


]} 


( 6 ) 


Note that eq. (6) can be rewritten in terms of galactic latitude and longitude (b,l) with 
the identification cos0 = cos 6 cos/. The next step is to determine what fraction of the 
halo’s mass is made up of axions. 

In clusters of galaxies I will argue that the axion mass fraction should be > fi a . This is 
not the case in the halo of our galaxy. In our halo, phase space arguments similar to those 
proposed by Tremaine and Gunn 24 play an important role in limiting the axion density. (I 
will discuss this point in detail in the next section; see also ref. 25) Inserting the central 
halo density implied by eq. (5) into eq. (12a), the axion mass fraction is seen to obey 

r a (max ) ~ 1.1 x 10 -6 g a a 2 1+ ~^ 


where the halo velocity dispersion is a — cr 2 x 100 km s l . r„(max) is much less than £l a , 
implying a much weaker signal than for clusters of galaxies. (As I will show later, the Bose- 
Einstein nature of axions will do little to change r„. Furthermore, since m a cc I a , an 
uncertainty in r 0 (max), which itself varies as m 4 , of even a factor of 10 will not significantly 
affect the conclusion.) Making the identification = r a (max)E and using expressions 
derived in section III leads to the conclusion that axions in our halo should produce an 
emission line at a wavelength of A tt = 24800 A( eV/m«), with a width of ~ 20 ka 2 { eV/m.), 
and with an intensity of 


jHalo 

a 


(R, A) = 1.7 x 10~ 25 




-(A - A a ) 2 
A 2 



( 8 ) 


in units of erg cm -2 arcsec -2 A -1 s -1 . Here I have defined G(x,9) = x 2 /(l + x 2 )J(x,0). 
The function G(x,9) varies between about 0.1 and 1 for reasonable choices of x (1/4 < 

x ^ i). Ha[o 

The most simple-minded limit which can be derived from eq. (8) is to require that I a ° 0 
not exceed the continuum level of the night sky (Ins ~ 10“ 17 erg cm 2 arcsec 2 A 1 s x ). 
Assuming x ~ 1 (G ~ 1), I find the very conservative bound m a < 6C -1/S eV from the 
fra In of our galaxy. Any axion more massive than this would create a line in the spectrum 
of the night sky which could not have escaped notice. Thus, as claimed, the signal which 
would arise from axions in our halo complements the bounds arising from the DEBRA 
measurements. 
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Astrophysics and cosmology have provided several convincing ways of constraining the 
axion mass to delineate the 3 eV to 8 eV axion window. The lower bound is firmly placed 
by considerations of the neutrino burst of SN 1987A. Stellar evolution arguments may 
have already closed the window, but this is far from certain. The relic decays of highly 
clustered and unclustered axions provide a complementary, and conservative, probe to the 
upper bound of the window. Both of these methods of constraining the axion mass share 
a dependence upon the axion model. But even in the worst case scenario, E/N = 2, relic 
decays still constrain the axion to be less massive than about 8 eV. This constitutes the 
claimed upper bound to the axion window. In the remainder of the paper I will describe 
an attempt to probe this window using the decays of axions residing in clusters of galaxies. 

III. Axions in Clusters. 

In order to search for relic axions in clusters of galaxies it is important to carefully 
determine how the axion signal will manifest itself. To do this one must determine the 
abundance of axions in the cluster, how they are distributed, the shape of the line, and 
several other points. In this section I will address these issues and derive expressions for 
the expected signal from axion decay in clusters of galaxies. 

The first step in determining the axion abundance in clusters of galaxies is to find 
their universal abundance. (This was important in deriving the previous bounds to the 
axion mass from the DEBRA flux. In this argument I adopt high energy physics units 
h = c = k B = 1.) As mentioned above, axions with masses greater than 10 2 eV are 
thermally produced. 10,12 The relevant process for creating a thermal population of axions 
is axion-pion conversion: N +ir — + N +a where N is a nucleon of mass — 939 MeV . The 
reaction rate for this process is roughly given by T = n/v(<r|u|) /■**/ (TVm’JlmjvA.WJV)] 2 
(m^/T) -3 / 2 exp(— m^r/T). The rate to compare this with is the expansion rate, H = 
1 .66gl^ 2 T 2 /m p i where m p i is the Planck mass and g» is the number of relativistic degrees 
of freedom. 26 When T/H > 1 for an expansion time, (~ 1 /H) axions thermalize. The 
ratio of the two rates is 

r/ ff ~io»(^) 2 (^)- 3/ %xp ( - mw /r). (9) 

For axions with masses of several eV, T/H is larger than 1 for temperatures above 50-60 
MeV. Thus at temperatures above 60 MeV axions are in thermal equilibrium. (Note: the 
above expression is only valid for temperatures, T < Tq/h ~ Aqcd — 150-400 MeV. 
This is the temperature of the quark-hadron phase transition and before this time there 
are no nucleons or pions to interact with. Axions will likely be in equilibrium above this 
temperature due to their interactions with quarks, but here we are only interested in their 
freeze out temperature which lies below Tq/h.) At temperatures below 50 MeV the axion 
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abundance freezes in at a value of n a = (C(3)/7r 2 )T 3 where T a is the axion temperature. 
This corresponds to a present day axion number density of 


n„ ~ 55 cm 



( 10 ) 


where g»F is the number of relativistic degrees of freedom when axions froze out at 
T ~ 50 MeV . At temperatures of ~ 50MeV the average energy per particle is ~ 
2T = 150 MeV , thus the particles which should still be relativistic are the 7r°, e , 
v v. ... 7 a. Naively summing these degrees of freedom results in <7*F = 18.75. This 
is not quite correct, since both the pions and the muons are becoming non-relativistic at 
this temperature. The result of a numerical integration 26 gives g. F ~ 15, the value I use 
throughout. This can easily be converted into a mass density and thence to a fractional 
density of the Universe 

( 11 ) 


^0^5° — 49 ( e v) L. f ) 


and Ho = 50/i 50 km/sec/Mpc. Eq. (6) shows that the axions in the multi-eV window 
cannot be the dark matter which seems to pervade the Universe, but that they do con- 
tribute an am ount of density comparable to that of baryons. Having determined it is 
now important to determine what fraction of galaxy cluster’s mass could be in axions. 

Axions, along with baryons and any other non-relativistic particles, should partic- 
ipate in the collapse which forms clusters of galaxies. The equivalence principle en- 
sures that the ratio of the axion number density to the number density of all species 
which collapse into the cluster should remain constant. This leads to the relation that 
Pa = (fi«/ ^cluster )pdust e r where p a is the mass density in axions, Cluster is the mass 
density of all matter in the cluster (which is determined observationally), and ft c luster is 
the fractional density of all matter that falls into the cluster. Since baryons are known to 
exist and should collapse, we know that ^cluster >&a + &B- An upper bound to ^cluster 
is ftciuster < ^TOT = 1- Thus the density of axions in the cluster lies between the values 
Pcluster >p a >n a Pcluster . Combining this argument with the well known bounds 
from primordial nucleosynthesis, 27 0.044 > QbH 0 > 0.15, we find that the axion density 
lies in the range 10 Q a p c i uste r Z Pa Pcluster for all axion masses in the multi-eV win- 

dow. I s hall use the most conservative assumption throughout the rest of the discussion, 
Pa = fta Pcluster , but it bears remembering that this may underestimate the number of 

axions in the cluster by up to a factor of 10. 

A possible loophole in the above discussion, which would tend to lower the number 
of axions in the cluster, is the question of available phase space. Once axions decouple 
from the plasma, their microscopic phase space density is conserved (i.e., they obey the 
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collisionless Boltzmann equation, df/dt = 0). For a thermal distribution of axions the 
initial phase space space occupancy is < 0(1) per cell of phase space for the great majority 
of axions. This is easily demonstrated: axions obey Bose-Einstein statistics. The axion 
distribution function is then f(E ) = {e E/T - l) -1 (ft = c = ks - 1), which is greater than 
one for E < E, = ln(2)T. The number density of axions with energies less than or equal 


to E. Is 


n(< E, 




E 2 dE 


e E/T _ ! 

where I have neglected m a (which is certainly valid at the decoupling temperature of 
~ 50 MeV ). For E/T < 1 the exponential in n(< E) can be expanded with the result 
that n{< E) ~ (g a T/2ir 2 ) EdE = (g a T /2it 2 )E 2 . This can be rewritten as n(< E) ~ 
(2C(3)) -1 (E/T) 2 n a ~ (2^(3)) -1 /(E) -2 n a , valid for E < T. We see that the number of 
axions which have phase space densities greater than 1 is suppressed by a factor of (E/T) 2 
at decoupling and that the larger the phase space occupancy, the smaller the number of 
axions with that occupancy. As an example of this, the above shows that the number 
density of axions with f(E) > 1 is ~ 0.2n o , or only about 20% (see also ref. 25). Carrying 
this even further, the percentage of axions with a phase space density greater than 10 per 
cell is only 0.3%. Furthermore, during the collapse of the cluster, or galaxy, dynamical 
processes such as phase mixing and violent relaxation are likely to reduce the phase density 
if they affect it at all. 28,29 

Since there is only a small fraction of the axions in high occupancy states and the 
phase space density is conserved, the arguments of Tremaine and Gunn 24 apply to thermal 
axions as well as to neutrinos. Assuming a phase space density of 1 and applying their 
argument to axions I find that the maximum axion mass fraction in a cluster based upon 
the available phase space, r a (moi), is 


r a (max ) 


(2t 


(12a) 


where a = <73 x 10 3 kms -1 is the one- dimensional cluster velocity dispersion and p c is the 
central density of the cluster. Going a step further and assuming that the cluster core is 
described by an isothermal sphere distribution or the analytic King model (to be described 
shortly) we can use the relation p c = §a 2 /±TtGa 2 to rewrite this as 


r a (max) 


2(2ir) 5 / 2 g a Ga 2 crm\ 

9ft 5 


~ 3 x 10 3 g a 



_ 2 

cr 3°250 


(12ft) 


where a = <2250 x 250 kpc is the cluster core radius. It is clear from the above discus- 
sion that the Bose character of the axion will not increase r a by very much. For axions 
with masses in the multi-eV window, with typical cluster parameters (<1250 = 03 = 1 ), 
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r„(max) > 1 f2 a . There is plenty of phase space for the axions in clusters of galaxies. 
This was not the case when we considered axions in the halo of our galaxy. In that case 
eq. (7) or a similar constraint plays a very important role in determining the axion signal. 

Thermal axions become non-relativistic when T a — m a /Z ~ C?(l) eV, before the time of 
matter-radiation decoupling (1 +z dee ~ 1150). Momentum, and hence velocity, are propor- 
tional to (1 f z )~ l f° r non-relativistic particles. Thus, unclustereu axions should be char- 
acterized by a velocity dispersion of order {v 2 ) 1 ^ 2 = v a ~ 4.3 x 10 4 c( eV/m„)(15/ <7*F ) ^ 
today. The requirement that axions collapse with the perturbations that form clusters (i.e., 
that the axion Jean’s mass is small enough) is simply that v e < d, the 3-dimensional 
cluster velocity dispersion, at the time of cluster formation. This requirement is amply 
met for cluster formation red shifts smaller than 15 or so. Thus, as claimed, there is every 
reason to expect that axions will collapse with other matter into clusters of galaxies. 

As the axions collapse, how will they distribute themselves in the cluster? The process 
of violent relaxation will tend to produce a Maxwellian velocity distribution. 28,29 ’ 30 Since 
this is the distribution function of an isothermal sphere one would expect the corresponding 
density profile. Numerical, N-body, experiments generally show that during collapse, a gas 
of dissipationless particles does not completely relax to the isothermal profile. The relaxed 
configurations are not too different from an isothermal sphere and this remains a useful 
approximation to the density profile. 31,32 (In fact, N-body simulations of cluster formation 
and relaxation result in profiles very similar to the analytic King, or modified Hubble, 
profile which I shall use to model clusters. 28,30,32 .) This picture is further bom out by 
observations of clusters themselves. Since clusters as a whole, as opposed to the individual 
galaxies within them, appear to have undergone little dissipation, the light should trace 
the mass in these objects; in any case, it should trace the gravitational potential well. It 
should then be possible to use both the galaxy and x-ray gas distributions to investigate 
the axion distribution. Observations reveal that the density distribution of the isothermal 
sphere and various approximations to it seem to provide an excellent description of the 
available data. 30 

A number of density profiles have been successfully fit to the distributions of galaxies 
and x-ray gas in clusters. 28,30 Among the most successful of these is the profile 


P( r ) 


Pc 

(l +r 2/ a 2)3/2 


(13) 


where p c = the central density of the cluster, r = the radial distance from the cluster 
center, and a. is the cluster core radius. This is often referred to as the analytic King (AK) 
model. The AK model has severed desirable features beyond its ability to fit the data. As 
its n am e implies, it gives simple analytic formulae for all of the relevant observables. It is 
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also a fairly accurate representation of an isothermal sphere for r < several a (the region 
where most of our data was acquired) and has a central density determined by observable 
quantities in the cluster 

Pc = 

where a is the 1-dimensional velocity dispersion of the cluster. Another profile which is 
often used to model clusters is 


9 a 2 


AizGa 2 


(14) 


p{r) = 


Pc 


(15) 


(l + r 2 /a 2 ) - 

This is a better fit to an isothermal sphere distribution in the outer regions of the cluster, 
but at the expense of the inner regions. When comparing the observations to the data, I 
have considered both of the above models. In this paper I will limit the discussion to the 
AK model, with only cursory comments about the model of eq. (15), for two reasons: 1) 
the AK model of eq. (8) seems to be slightly more successful in modeling clusters than the 
other distribution and 2) in the region which the observations cover, there is not a large 
difference in the predicted axion decay intensity between the two models; I obtain more 
conservative bounds using the AK model. Thus, any statement made concerning the axion 
based upon the AK model will only be stronger in the model of eq. (10). 

Having decided upon a density profile to describe the cluster, and hence the axion 
distribution in the cluster, the last ingredient needed to find the axion signal in a cluster 
is the line shape. Particles in the cluster are generally taken to be distributed with a 
Maxwellian (gaussian) velocity distribution 30 


V(v r )dv r = 


ry/^K 


exp[— (v r - (v r )) 2 /2(T 2 ]dv r 


(16) 


where v r is the measured velocity projected along the line of sight (recessional) of the 
particle (i.e., galaxy) in the cluster. Here V(v r )dv r is the probability that an individual 
particle has a lin e of sight velocity in the range v r to v r + dv r . This is easily converted into 
a probability distribution for A using the non-relativistic Doppler formula, AA/A = u/c, 
with the result 


V(\)d\ = 


■exp 


(A - A„) 2 c 2 


A 2 


2<r 2 


d A. 


(17) 


y/2ncr\ a 

Eq. (17) (without the d\) is the profile of the line. For future reference, eq. (17) describes 
a linp with a gaussian profile and a full-width at half-maximum (FWHM) of AA)fwhm = 
2\/ln4(<7/c)A a ~ 2.4(cr/c)A a . 

It is now straightforward to predict the expected intensity of an axion decay line as 
a function of A, projected cluster position, and axion mass. The first thing to do is to 
integrate along the line of sight through the cluster to obtain the surface density E(iZ). 
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(18) 


Using the AK model the surface density is 


X(R) = 


2p c a 

(1 + R 2 /a 2 ) 


where R is the projected radial distance from the cluster center. As per the discussion 
above, the axion surface density is E a (R) = Q a E(R). With this identification the predicted 
axion signal becomes 


t t r> \\ — = ^q(^) c2 ^( A ) 

M ' - dAdndXdt 47rTa 


(19a) 


Combining equations (2), (11), (14), (17), and (18) and putting in numbers this becomes 


'•(**) - 6 9 x 10 ~ 2 ' Gv ) 7 (ssd © c2exp ( 


'-(A-A a ) 2 c 2 


in units of erg cm -2 arcsec -2 A -1 s -1 . Eq. (19) describes an axion decay line, in the 
cluster rest frame, which is centered at a wavelength A a = 24800 A( eV/m a ) and a FWHM 
of ~ 200 Acr 3 ( eV/m a ). In our rest frame the line center is shifted to A a => A a (l + z) and 
the width is broadened by a factor of 1 + 2 as well. These two considerations will help to 
distinguish candidate axion lines from night sky lines. Finally, in our frame, I a suffers from 
a cosmological dimming which further reduces the flux by a factor of (1 + z) 4 . Important 
features to note about eq. (19b) are: the spatial dependence of the line strength (<x R 2 
for R > a), that it depends upon ( 2 , and the strong dependence on m a (to the seventh 
power). The dependence of I a oc p a /Kr a upon m a is easily understood, m\ from r a , m a 
from p a oc fi a , and m a from A a . Had we used the density profile of eq. (15), the main 
difference would have been in the spatial dependence, I a oc 1 / R for R > a. 

Now that we have a prediction for the axion signal, we can design a search which will 
optimize the chance of detecting this feature. Neglecting the line shape, and much of the 
parameter dependence, eq. (19b) can be rewritten as ~ 1.5 x 10 -17 (m a /3eV) 7 C 2 /(l + 
R 2 /a 2 ) erg cm -2 arcsec -2 A -1 s -1 . Thus an axion in the lower end of the mass window 
has a predicted intensity comparable to that of the night sky continuum (I N s ~ 10_17 
erg cm -2 arcsec -2 A -1 s -1 ). Furthermore, there are numerous strong emission lines 
throughout the relevant spectral range (~ 3000 A-8500 A ). 14 This problem is especially 
acute towards the red end of the spectrum (corresponding to small axion masses and 
hence to small I a ) where there are numerous emission bands from molecular OH (see Figs. 
2-6). In order to secure an unambiguous detection of an axion line, especially in the case 
of Q < 1, the night sky background must be removed. Fortunately, the spatial dependence 
of eq. (19) suggests an easy solution. By talcing a spectrum from near the cluster center, 
where the line is strongest, and subtracting from it a spectrum taken from the outskirts of 
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the cluster ( R > 2a), the night sky may be removed with very little effect upon the axion 
line. In principle, this procedure should allow the entire window to be easily probed. In 
practice, there are several complications. Foremost among these is the fact that the night 
sky varies both spatially and temporally, making a perfect subtraction impossible. In order 
to achieve the best subtraction possible one should to take the spectra as close together in 
space and time as possible. When this is done, intensities of between 1-10% of the night 
sky continuum are typically achieved. Another worry is that the axion distribution may 
be more diffuse than the baryon distribution (e.g. a ax i on >• a yo / ariea ) and the subtraction 
procedure may remove the line. I will address this concern in the next section. 

A search for an axion decay line in three well-studied clusters of galaxies A1413, A2218, 
and A2256 was carried out in May 1990. 15 These three clusters have a number of properties 
which make them ideal for a search of this type. Each of these clusters has been extensively 
studied at both optical and x-ray wavelengths. Thus, the cluster parameters such as the 
core radius, a, and the velocity dispersion, <r, axe fairly well determined. They axe quite rich 
clusters with well determined centers; each has a fairly large velocity dispersion, meaning 
a laxge mass and therefore a potentially stronger axion signal. Another advantage of these' 
clusters is their distance from us. Each is distant enough so that the cluster is reasonably 
small on the sky, 0(5 arcmin), yet not so fax away that candidate axion lines axe shifted 
too far into the red or dimm ed by too great a factor. Additionally, they span enough region 
of red shift space so that if the axion line falls in a region of poor night sky subtraction in 
one cluster’s spectrum, it should fall in a region of good subtraction in another’s spectrum. 
Table 1 summarizes the relevant parameters of the three clusters. There is some scatter in 
the data, so in choosing the values of a and a to use in eq. ( 19) I have chosen representative 
and convenient values which agree with most of the data. For the actual data see references 
33, 34, and 35 for the clusters A 1413, A2218, and A2256 respectively. 

Before proceeding to describe the observations of these clusters I will discuss a final 
issue which bears upon the multi-eV axion window. Now that we have a firm prediction 
of the intensity of an axion decay line in a cluster of galaxies, we can ask if any previous 
searches for intracluster light might have detected it. The answer turns out to be “yes” 
for a range of mass both within and outside of the window. In the ultra-violet region 
there are two spectroscopic observations of the Coma cluster of galaxies which clearly 
rule out axions in the mass ranges 16 eV-20 eV and 22.1 eV-27.8 eV for any plausible 
value of £. 36 As support for this statement I consider the data of Henry and Feldman 36 
at 1300 A ( m a = 19.5 eV). Their upper limit to any line-flux in the Coma cluster at this 
wavelength is 1.4 x 10 -19 erg cm -2 arcsec -2 A -1 s -1 . For an axion decay line with 

£ = 1 eq. (19b) predicts I a ~ 1.7 x 10 -12 erg cm -2 arcsec -2 A -1 s -1 ; a flux of more than 
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10 million times that measured. This strongly demonstrates the strong dependence of I a 
* upon the axion mass. To hide such a line in Coma, one needs C < 0.0003, requiring a 
near miraculous cancellation to the two amplitudes which lead to axion decay. Of more 
interest are two sets of broad band searches for intracluster light performed in the optical. 37 
Shipman and Cowsik 37 used a number of pre-existing observations of mtergalactic light 
to set constraints upon the lifetime of the neutrino to radiative decay. Their limi»o may 
easily be translated to limits on m a via a modified eq. (19). Since the observations are 
broad band, they only measure the total flux in the band; all wavelength information is 
lost. It is simple to predict the total flux produced by an axion line in a band; simply 
integrate eq. (19) over the wavelength, A. This has the effect of making the function 
V(\) become 1 with the corresponding change in I a . When this is done, the g and r band 
observations used by Shipman and Cowsik rule out the standard (C = 1) axion in the range 
5 4 e V-3.6eV. The work of Partridge and collaborators 37 is similar but probes further into 
the red. Unfortunately, due to problems with the data it is uncertain how to convert their 
data into constraints on m 0 . Previous data place interesting limits upon m a but do not 
close the axion window. FVirthermore, these limits are difficult to turn into convincing 
constraints on this is much simpler for spectroscopic observations, such as those I will 
describe shortly. 

IV. The Observations. 

Long slit spectroscopic observations of the clusters A 1413, A2218, and A2256 were 
obtained on the nights of 24 and 25 May 1990 at KPNO. The 2.1 meter telescope was used 
with the Gold camera CCD spectrograph at a resolution of ~ 10 A. This resolution is well 
matched to a cluster axion line which would have a FWHM of about 200<73 A(eV/m a ) in 
the three clusters. The incoming light passes through a 5' by 2."5 slit onto the grating 
of the spectrograph which then disperses the light onto 480 x 800 pixels of the CCD. 
Each exposure contains both spatial (480 pixels at 0. M 78/pixel) and spectral (800 pixels at 
4.6 A/pixel) information. Pixels were binned by three in the spatial direction to increase 
the signal to read noise ratio in the CCD. Thus, there are 160 individual spectra per slit 
position, each sampling a different spatial point in the cluster. On the first night we took 
multiple exposures of all three clusters using KPNO Gold camera grating #400. This 
grating allowed observations in wavelength range 4762-8441 A, corresponding to axion 
masses in the range 6.1-3.2eV when red shift effects are taken into account. For each 
cluster we w alke d the slit from a position near the cluster core to the cluster outskirts 
along an E-W axis, obtaining spatial information out to between 3 and 10 core radii 
for each cluster. Thus, if there were an axion line present, its spatial profile could be 
determined. This would help to confirm the line as arising from axion decay as well as 
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providing valuable information about the cluster potential. On the second night grating 
#201 was used, which allowed a wavelength coverage of 3737-7606 A (axion masses from 
7.8-3.7eV, including red shift), to observe the two smaller clusters (A1413 and A2218). 
Here single exposures were taken of each cluster as the 5' slit covered more than 3 core 
radii for each. The exposure times ranged from 30 to 75 minutes per slit position and were 
chosen to ensure a large signal „o-noise ratio for the postulated axion line (see table 2 for 
observing information). 

The data were reduced using the Image Reduction and Analysis Facility (IRAF) version 
2.8 as supplied by KPNO. The CCD bias was subtracted from each frame, and spatial 
variation in the pixel response and illumination were removed by dividing each spectrum 
by dome and sky flats taken earlier in the evening. Calibration in the spectral dimension 
was obtained by exposures of a He-Ne-Ar lamp. Observations of a KPNO standard star 38 
(Wolf 1346) provided absolute flux calibrations for all exposures. 

In Figs. 2-6 I present the spectral data used in the search for an axion decay line. 
Figs. 2-4 are from the first night’s observations and 5 and 6 are from the second night. 
In each of these sets of data (a) is the inner aperture, or “on-cluster” data. The on- 
cluster data is the spatial average of 30 pixels (23. "4) from near the cluster core, R/a < 1. 
Part (b) of each figure is the “off-cluster” data. The off-cluster spectra are the average 
of 30 pixels well away from the cluster center, R/a > 3. Comparing Figs. 2-6 it is 
immediately apparent that all of the spectra are basically similar; they are simply spectra of 
the night sky. General features worth noting are the relatively constant continuum level of 
Ins ~ 10 -17 erg cm -2 arcsec -2 A -1 s -1 and the numerous emission lines and bands which 
stand high above the continuum level and vary in intensity from spectrum to spectrum. No 
line which meets the criteria for an axion decay line, set out below, is in evidence. To lessen 
the chance of inadvertently removing a candidate axion line, cosmic-ray hits on the CCD 
have not been removed. These hits are easily distinguished on an individual basis as they 
meet none of the criteria which define an axion line. The bright level of the continuum and 
the many lines could very well overwhelm, or mimic, a candidate cluster emission line, so 
it is advantageous to remove them as well as possible. This is accomplished by subtracting 
the off-cluster spectrum of (b) from the on-cluster spectrum of (a). The resulting “on-off” 
spectra are shown in Fig. 7 a-e. This procedure should efficiently remove the night sky 
and, by virtue of the spatial dependence of the axion line (see eq. (19)), scarcely affect 
any cluster emission fine at all. 

The most striking fact about the on-off spectra in Fig. 7 is the quality of the subtrac- 
tion. Intensity levels of 1-10% of the continuum night sky level are easily achieved, except 
in the very brightest of the night sky emission lines. Furthermore, there appears to be 
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little or no offset of the subtracted spectra in Fig. 7 from the value of 0, again indicating 
the quality of the subtraction. The best subtraction achieved is that in Fig. 7b, A2218 
on night 1. Because of the relative smallness of the cluster, a = 0/88, compared to the 
size of the slit, ~ 5', both the on-cluster and the off-cluster spectra came from the same 
observation, eliminating the temporal variation of the night sky. This spectrum has several 
features which occur in all of the on off spectra and so deserves a bit of discussion. While 
in general the subtractions seem to be quite successful, there appear to be a number of 
residual lines — might these be the result of axion decay? For this spectrum, and all of the 
others, the answer is “no”; all of the residual lines fail to meet the axion line criteria. In 
Fig. 7b the two lines, one positive, one negative, at A = 5035 and 5350 A are cosmic-ray 
hits. The features at 5577 A and to the red of 7500 A are due to poor night sky subtraction 
in the bright emission lines (and bands). These features are endemic of all of the spectra 
and each must be checked individually. One last point to note about the spectra in Fig. 
7 is the relatively large amount of noise in the blue portion of the second night s data. 
This is due to too little flux reaching the pixels in our dome flats in the blue end of the 
spectrum (the spectrograph’s efficiency drops rapidly for short A). This noise turns out 
not to be a problem because of the strong dependence of I a on m a (/ a oc m a ). Any line in 
this part of the spectrum would stand far above the noise for any reasonable value of (. 

What criteria must a candidate line in the spectra of Figs. 2-7 meet before it can be 
seriously considered as being from axion decay? First, each of the clusters should show a 
line that has been red shifted from a common rest wavelength. Next, the line should have 
approximately the predicted intensity, eq. (19), and be Doppler broadened to a gaussian 
with a FWHM of AA/A = 2\/ln4<x/c. Finally, the line should have a reasonable spatial 
profile (i.e. follow the AK or some similar model). Besides using the data in Figs. 2-7, this 
could be checked by examining spectra from apertures between the two already examined. 
Not only would this be a useful check on whether the candidate line could be from axion 
decay, but could yield interesting information on the cluster potential. While these criteria 
are not enough to make the case for an axion decay line with certainty (eg. an intracluster 
atomic emission line could easily satisfy all of the requirements), any line meeting them 
would certainly bear closer scrutiny by other methods, (e.g. a proposed axion-photon 
conversion experiment which could detect solar axions in the range 0.1 eV ^ tn a iS 5 eV. ) 
Naturally the converse is true; any line not meeting these requirements can be discarded 
from consideration. 

The simplest way to check for the presence of an axion line in the data is to compare 
the real data of Fig. 7 to data with an artificial line, which obeys eq. (19), inserted into it. 
In Fig. 8 I show several spectra with axion decay lines which arise from axions of different 
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mass. Fig. 8a, when compared with Fig. 7c, illustrates the major result of this paper — that 
there is no intracluster line emission from axion decay for axions with m a > 3.2 eV. Fig. 
8a is the on-off (R/a = 0.48 minus R/a = 2.96) spectrum of A2256 with a 3.2 eV axion 
line artificially introduced. The line is at a central wavelength of A = A a (l z) = 8216A 
with a FWHM of ~ 80 A and an intensity given by eq. (19b) (with ( = 1). The axion 
line is the most obvious feature in the ..pectrum and is easily identifiable, in spite of the 
fact that it lies in the noisy, red portion of the spectrum. Since /„ oc m 7 , any line from 
a more massive axion would stand out far more prominently. A2256 has the lowest red 
shift of the three clusters surveyed and hence probes the lowest range of axion mass; Fig. 
8a represents the lower bound to the axion mass from this search: m a < 3.2 eV. The 
remainder of the data will serve to confirm this result, as well as exclude the presence of 
an axion with m a < 7.8 eV. In Fig. 8b I have placed a 3.5 eV axion decay line in the night 
1 spectrum of A2218, and in Fig. 8c a 6.5 eV line in the night 2 spectrum of A2218. This 
latter figure should be compared to the spectra of Fig. 6; a line this intense would stand far 
above the continuum in the unsubtracted spectra as well. The important point to notice 
in these figures is the strength of the line compared to the zero level of the subtraction. 
Such lines axe clearly not present in any of the spectra of Fig. 7. The final result of these 
comparisons is that an axion with a mass in the range 3.2 eV < m a < 7.8 eV is firmly 
excluded by the observations, effectively closing this window of axion mass (for £ = 1). 

Having now established that the existence of the “simplest” axion ( E/N = 8/3, £ = 1) 
is excluded by the data, what additional information can be extracted from our data about 
axions arising in more complicated models (£ < 1)? The model dependence will turn out 
to be slight, since any bound to the axion mass based on eq. (19) will vary as m a <x £ -2 / 7 . 
Since no cluster line is in evidence in Fig. 7, we cannot simply read off the value of £ 
by comparing the line intensity with that predicted. The best alternative to this is to 
determine the limit which the data places on any line emission in the cluster and, via eq. 
(19), use these limits to determine the minimum value of £ that would produce a line at 
or below this intensity level. This problem will be attacked in two manners. First I will 
discuss a cross-correlation technique where the spectra of Fig. 7 are correlated with each 
other and noiseless “template” spectra in search of a correlation peak which would indicate 
cluster line emission. The cross-correlation technique is well suited to detecting a cluster 
emission line but is too time consuming to establish absolute line-flux limits throughout 
the spectrum. The second method I will discuss is simply to ran a moving window, with 
a width equivalent to an axion decay line, across each spectrum and calculate the mean 
flux in each window. The “two-sigma” upper bound on this flux is then a firm bound on 
any intergalactic line-flux in the cluster and can be used to constrain £. 
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The use of Fourier cross-correlation techniques to pull signals out of the noise is often 
used in many spectroscopic applications. 40,41 In searching for intracluster line emission the 
application of this technique is simple and straightforward. While the wavelength of the 
would-be line is unk nown, a line intrinsic to the cluster will appear at different wavelengths 
in different clusters, Aj = A a (l -h *,). Upon taking the logarithm of this equation we find 
In A; = In A a + ln(l + Z{). If the spectral dimension of the spectrum is converted to the 
logarithm of A, lines in different clusters will be offset by a constant, known amount no 
matter what the intrinsic wavelength of the line. Since this is the case, an intracluster line 
will lead to a positive peak in the cross-correlation function of the on-off spectra for any 
two clusters 

_ f J 1 (x)f 2 (x + /)di _ f h(k) m l2(k)acp(-ikl)dk ^ 

{ ~ (fifjx jqdx) 1 ' 2 (/|/i(t)l 2 * /|i2(*)l 2 *) 1/2 ’ 

at a lag / = ln[(l + z 2 )/(l+z l )}, where the variable is x = ln(A). Moreover, the height and 
width of the peak in £(/) are directly related to the intensity and width of the intracluster 
emission line. A statistically significant peak at the correct lag and of the proper width 
would provide strong evidence for an intracluster line. If such a peak was found in «0 the 
spectra could then be carefully searched in order to determine A,-, the wavelength of the 
line, and then m a . 

To assess the statistical significance of a peak, one must know the distribution of false 
peaks due to noise. In each of our on-off spectra all of the prominent features which are 
obviously associated with cosmic-ray hits or poor night sky subtraction have been removed. 
In Fig. 9 I show the on-off spectrum of A2256 with major cosmic-ray hits and night sky 
lines removed (compare this with Fig. 7c). The central limit theorem then ensures that 
the re mainin g noise peaks have a gaussian distribution. Hence the probability of a positive 
noise peak of height greater than h is 

P{> h) = j exp(-h?/4crt)dh/y/TT(Ts = eric(h/2as), (20) 

where is the rms of the antisymmetric part of £(/). In computing <r $ one must eliminate 
the possible contribution of real signals. This is done by calculating the rms value of the 
asymmetric part, £asym(0 = + /) — 0 } o — Ol/^i where Z<j is the lag expected for a real 

signal, and ass umin g that the rms value of the symmetric part is the same. Any real signal 
doesn’t contribute because its asymmetric part vanishes. A more detailed discussion of 
this technique and the underlying theory is given in reference 40. 

In Fig. 10 a & b I show 2 cross-correlation functions between on-off spectra of clusters. 
Fig. 10a is the night 1 spectra of A2218 and A2256 cross-correlated; Fig. 10b is the night 
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2 spectra of A1413 and A2218 cross- correlated. For comparison, in Fig. 11 I include two 
correlation functions with artificial axion lines introduced. In Fig. 11a I correlate the night 
1 on-off spectra of A2256 and A2218 with a 3.5 eV axion decay line added to each. In 
Fig. lib I include a 6.0 eV axion decay line with ( = 0.07 ( E/N = 2) in the night 2 on-off 
spectra of A2218 and A1413. In both cases the correlation peak clearly stands out. 

The cross-correlation is carried out using the rvxcor task in IRAF. This task takes 
the two spectra to be correlated and returns £(/) (as in Figs. 10 and 11), the lag of 
the correlation peak in pixels (which can then be converted into A), the normalized peak 
height h (0 < h < 1), and the value of cr^. Using data with artificial axion lines introduced, 
it is simple to determine the correct lag of the peak (in pixels) for the correlation of cluster 
with cluster. Having determined the location of the presumed peak, a gaussian is fit to this 
region of £(l) for the unaltered spectra to determine both h and cr^. These two quantities 
can then be used to determine the statistical significance of any peak present. The criteria 
which a “real” correlation peak must satisfy are that it must be within l-<r of the correct 
lag and it must have less than a 5% chance of being a noise peak (V(> h ) < 0.05). This 
corresponds to a 2 -a or greater probability that the peak is due to a real correlation and 
is not a spurious noise peak. No peak meeting these criteria is present in the correlation 
functions of Fig. 10. As an example, when a fit is done to the correlation peak nearest 
the correct lag for the night 1 data of A2218 and A2256, the fitted peak is at a lag 15% 
too small and has a height in the noise-peak distribution where there is a 39% probability 
of having a peak of this height or larger. This peak is clearly not from an intracluster 
emission line. 

To set flux limits using the cross-correlation technique, I have modified the procedure 
somewhat. Instead of correlating the cluster spectra with one another, I correlate one 
of the cluster spectra with an axion decay fine inserted against a noiseless “template” 
spectrum containing the same mass line (the template is simply the fine described by eq. 
(19b) with nothing else added). The strength of the line is then reduced until the height 
of the cross-correlation peak is such that there is only a 5% chance that the peak is due to 
the gaussian noise. The intensity at which this happens then constitutes the 2 -a limit to 
any line-flux at that wavelength. These limits on the flux may then be turned into limits 
on the model parameter £. This procedure has been carried out for a number of axion 
masses using the data; the most restrictive limits are summarized in table 3 and Fig. 13 
by the dashed line (these results were also presented in ref. 15). 

The cross-correlation method discussed above produces quite stringent limits to the 
flux of intra-cluster light and hence on the parameter (. An examination of Fig. 12 
reveals that axions more massive than 4 eV seem to be excluded by the data, even for 
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E/N = 2 (( = 0.07). The disadvantage of this technique is that finding absolute limits for 
every value of m a is time consuming; one must construct template spectra at the central 
wavelength of the 796 pixels of each spectrum and cross correlate them with the spectrum. 
As an alternative to this, I now investigate another way of getting limits to any line-flux 
in the on-off spectra. 

The spectra of Figs. 7a-e represent useful limits to the presence of any line-flux in the 
three clusters of galaxies observed. One method of quantifying the amount of flux in each of 
the spectra is to sequentially examine a series of small “windows” in the spectra, each with 
a width roughly equivalent to that of an intracluster emission line. The FWHM of an axion 
line intrinsic to a cluster of galaxies is 210 A<T3( eV/m a )(l + z). For axions in the multi-eV 
window this line width will vary between about 40 A and 90 A. This suggests that by 
examing the flux in a small window of the spectrum of width 90 A (20 pixels) a limit to the 
flux of an axion line in this window can be obtained. In each window, the mean intensity 
per pixel, 7, and the standard deviation from the mean per pixel, 07, can be calculated. 
The total integrated flux in the window is then given by Fw = (4.6 A /pixel) x 20pixels x I 
with a standard deviation given by a w = (4.6 A/pixel) x v^Opixels x 07. a w is used to 
determine the distribution of noise in the spectrum. We can now use Fw and a w to set 
quantitative limits to any line-flux in the window by comparing Fw with the integrated 
flux that a cluster emission line would produce in the window. The probability that a 
cluster emission line with a flux, F a , would produce a detected flux of less than Fw in the 

datais 1 rFw —Fa 

V{F a ,F w ) = -J=^ / exp(—h 2 /2cr 2 )dh, 

A/27T0- 7-00 

where the distribution of noise peaks is taken to be gaussian. (This expression derives 
from requiring that the noise, when added to the signal (F a ) results in a detected flux 
of of Fw or less.) In order to set 95% confidence limits to the presence of any line-flux 
(similar to those derived via cross-correlation) I require that only 5% of the time is a 
downward fluctualtion of the noise sufficient to suppress a signal to a flux less than that 
observed, F w (i.e. V{F a ,F w ) = 0.05). This occurs when F w - F a = -1.65 aw- Thus 
any cluster emission line with a flux 1.65&W greater than that observed in the window 
can be excluded at the 95% confidence level. Now that we can calculate a limit to any 
flux in a 20 pixel window we can move such a window throughout an entire spectrum, one 
pixel at a time, and derive limits to the amount of line emission in each of the subtracted 
spectra of Figs. 7a-e. In Figs. 12a-e I present the 95% confidence limits to any line-flux 
in the observed clusters by this technique. In Figs. 12d and e, the large downward (and 
negative) fluctuations in the blue are due to the problems with the dome flats alluded 
to earlier. Note that the flux levels in Fig. 12 are comparable to those obtained via the 
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cross-correlation analysis. 

The flux limits displayed in Fig. 12 provide a means of determining firm limits to £, 
the axion model dependence parameter. By comparing each of the 5 on-off spectra at 
points corresponding to a single axion mass and by choosing the lowest value of the flux 
at this point, the maximum value of /„ consistent with the observations can be identified. 
F-. (19) can then be used to determine the maximum value of £ which produces an axion 
decay lin e with an intensity equal to the observational limit. The results of this exercise 
are displayed in Fig. 13. The two bumps in £ near m a = 7eV arise from the overlap 
of the regions of negative flux in Figs. 12d and 12e. Since the data is unreliable here 
I have simply taken the flux limit to be the level of the night sky continuum, Ins = 
10~ 17 erg cm - 2 arcsec - 2 A - 1 s -1 . The values of £ are still less than for the E/N = 2 
case. It is clear from this figure that the maximum axion masses allowed for E/N = 2 
are m„ < 4eV; far below the equivalent bounds from RG-HB stellar evolution. Finally, 
it is worth remembering that I have chosen the mass fraction of axions in clusters very 
conservatively (= Q a ). The mass fraction could be up to a factor of 10 larger. If this is the 
case, the limits to £ would improve by a factor of about 3 since, for a given m a , £ oc \fT~ a . 

As a final topic, I will discuss the possibility that the axions are not distributed like the 
luminous matter in the cluster. All of the limits considered so far assume that the axions 
have a spatial distribution similar to the galaxies and x-ray gas in the cluster, £ a oc 1/R 2 
or 1 /R for R > a, where a is determined by the luminous material. While this seems, by 
far, the most reasonable assumption, it may not be the case. Axions could, for example, 
have a similar spatial distribution but with a ax ion ^ a galaxies - This could occur if, for 
example, the majority of the mass in the cluster, M, consisted of cold dark matter (CDM). 
West and Richstone 42 have shown that, for Me dm > M ga [ ax i e3 , dynamical friction may 
lead to a strong segregation of the galaxies to the central regions of the cluster (this 
could also affect the x-ray gas if it is “blown out” of the galaxies after cluster formation). 
Since the axions would be distributed more diffusely than the baryons, (in one model of 
ref. 42, a ax i on3 ~ clcdm ~ (5-10)a J?a / a iie3) the predicted axion signal would be smaller 
( I a oc 1 la ax ion), and the subtraction would almost totally remove the line, significantly 
altering the limits quoted. 

This concern can be addressed by subtracting the data of one cluster from that of 
another cluster. Since any line due to the sapae physical process in the two clusters will 
be red shifted to different wavelengths, the line will not be removed by the subtraction. 
In fact, this procedure produces a unique signature for a cluster emission line — a pair of 
lines, one positive and one negative, corresponding to a single rest wavelength. In Fig. 14 
a h b I show two such spectra: a) the inner aperture of A2256 minus the outer aperture 
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of A1413 and b) the inner aperture of A2218 minus the outer aperture of A2256. Both are 
from the first night’s data. These two sets of spectra were chosen because they were taken 
close together in time, if not in space. There are two important features to notice in the 
on-off spectra of Fig. 14. First, the night sky subtraction is significantly poorer than the 
same cluster subtractions due to the large spatial separation of the apertures (a factor of 
~ 4 wcise, as quantified by the moving window technique). Note also, that the baseline 
value of the subtractions is once again close to zero (although not as convincingly as in 
Figs. 7a— e). Second, there is absolutely no indication of the existence of any intra-cluster 
line emission in either of the spectra. Since the axion core radius is being treated as an 
unknown here, no real predictions of I a can be made for the clusters. Hence, no statement 
concerning limits to £ can be made based upon the data if the axions do not track the 
light. 

V. Concluding Remarks. 

Astrophysics and cosmology have much to say about the axion. In this paper I have 
considered the relic decays of axions to photons in a number of astrophysical settings. 
By considering the effects of multi-eV axions diffusely spread throughout the cosmos and 
tightly clustered in the galactic halo, an upper limit to the axion mass of ~ 8 eV is obtained, 
even with the most pessimistic model dependence ( E/N = 2). In concert with the axion 
mass bound from SN 1987A ( m a > 3eV), these considerations lead to a well defined 
window of allowed axion mass, 3 eV < m a <8 eV . 

This window has been unsuccessfully searched by looking for an emission line from 
axion decay in the intergalactic light of three clusters of galaxies. No signal matching that 
expected from axion decay is seen in the clusters, ruling out the existence of the axion (in 
the simplest models, E/N = 8/3) in the mass range 3.2 eV to 7.8 eV. This implies, that if 
the axion exists, its mass lies in the range 10 3 eV — 10 6 eV and that it likely comprises 
the majority of the dark matter known to exist in the Universe. Limits to the flux of any 
intracluster line emission were obtained by two techniques: by cross-correlating spectra 
and by measuring the mean flux, in a moving window throughout the spectra. It should be 
stressed that these line-flux limits are very dependent upon knowing the spatial distribution 
of the emitting source, axions or otherwise. It is highly likely that the intergalactic matter 
should track the luminous matter (galaxies and x-ray gas) and this assumption has been 
used throughout. Under this assumption, the flux limits lead to interesting limits upon the 
axion model dependence parameter (j (or E/N). By insisting that the flux of any axion line 
be less than the observed flux in the cluster it is found that even for C = 0.07 (E/N = 2), 
the “worst case scenario,” axions with masses greater than 4 eV are excluded by the data. 
Finally, the question of whether the galaxies and gas in the cluster trace the mass has been 
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addressed by subtracting one cluster’s spectrum from another’s spectrum; no evidence for 
line emission is seen. The remaining astrophysical uncertainty of note is the mass fraction 
of axions in the clusters. This fraction has been treated very conservatively and could be 
up to a factor 10 larger than has been assumed. If the axion mass fraction in clusters 
was ~ 10f2 a , instead of fl a , axion masses significantly below 3eV could be probed by this 
method (a.-. axion of m a = 2.3 eV would produce a decay line with the same intensity as 
the 3.2 eV line predicted by eq. (19b)). 

A last concern I will address is the possibility that the window is not quite closed. I 
have used throughout the lower bound to the window of 3 eV. I have also stated that the 
data only extends to axion masses of 3.2 eV. (The data actually extends to m a = 3.1eV 
but the noise is too large to convincingly rule out axions this massive, although there is 
no evidence for one.) Could the axion lie in this last small bit of parameter space? This 
concern could be exacerbated if the uncertainties in the SN 1987A analysis all tended to 
push that limit in the direction of lower axion mass. The answer to the above question is 
“yes, it might” although the remaining parameter space is quite small. Unfortunately, due 
to the large amount of noise in the red portion of the spectrum this question is unlikely to 
be answerable by ground based observations. Another method would be needed to probe 
this range of axion mass (e.g. the proposed experiment of ref. 39), or the Hubble Space 
Telescope (HST) might answer the question once it is fitted with instruments which see 
farther into the red. High above the Earth’s atmosphere HST would be free of the variable 
night sky and could thus obtain much better subtractions, allowing the region m a < 3.2 eV 
to be probed. As an example, if HST’s subtractions red-ward of 8400 A were comparable 
to the best subtractions obtained in the blue portion of this investigation, line intensities 
of ~ 10 -18 erg cm -2 arcsec -2 A -1 s -1 should be detectable. (This value is typical of the 
limits obtained by the cross-correlation analysis in the quiet portions of the spectrum, see 
table 3.) One could then probe down to axion masses of ~ 2.5eV(<C = 1), corresponding 
to A = 10600 A in A2256. 

As a final comment upon this work I mention that, while everything done here has been 
specific to the axion, all of the results can be carried over to any relic particle which decays 
to optical photons. Equation 19 can be generalized to another relic, X, by multiplying it 
by (n- r /2)(£lx/tta)(T a /Tx), where n 7 = 1 or 2 is the number of photons produced in each 
A” decay. The limits which this data places upon the lifetime of a relic neutrino which 
decays into a lighter neutrino and a photon will be considered elsewhere. 15 
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Table 1. Summary of properties of the three clusters observed. 


Cluster 

a 

(km s -1 ) 

a 

kpc (arcmin) 

z 

A1413 33 

A2218 34 

A2256 35 

1230 

1300 

1300 

400 (2.03) 

200 h^j 1 (0.88) 

473 h" 1 (5.0) 

0.1427 

0.171 

0.0601 


Table 2. Summary of information concerning the observations of the three clusters ob- 
served: exposure time is per slit position, if more than one slit position was observed for 
a cluster the slits lie on an East- West axis with the slit centers separated by 253", night 1 
refers to 24 May and night 2 to 25 May. 


Cluster 

Mass 

range 

(eV) 

Exposure 

time 

(sec) 

Inner 

aperture 

(fi/a) 

Outer 

aperture 

(fi/a) 

Number 

of 

Slits 

Night 

A1413 

3.36-5.95 


1.11 

4.64 

3 

1 


3.73-7.58 


0.65 

2.94 

1 

2 

A2218 

3.44-6.1 


0.94 

5.33 

2 

1 


3.82-7.77 


0.94 

5.33 

1 

2 

A2256 

3.11-5.52 


0.484 

2.96 

3 

1 
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Table 3. Summary of limits to the flux and the parameter C from the cross-correlation 
technique. 


m a 

Flux 

c 

Cluster 

3.5 

1.35 x i0" 18 

0.16 

A2256 

4.0 

8.0 x 10" 19 

0.078 

A2218 

4.5 

4.6 x 10“ 19 

0.039 

A2218 

5.0 

6.6 x 10“ 19 

0.032 

A2218 

6.0 

5.9 x 10" 19 

0.016 

A2218 

7.5 

1.3 x 10 -18 

0.011 

A2218 
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Figure Captions 


Figure 1 Constraints to the parameter ( obtained by considering the effects of the electromag- 
netic decay of unclustered axions upon the DEBRA. 

Figure 2 Spectra of two intergalactic regions in the cluster A1413, night 1. a) The inner, or 
“on,” aperture (R/a = 1-11); b) the outer, or “off,” aperture ( R/a — 4.64). Intensity 
is in units of 10 -18 erg cm -2 axcsec -2 A -1 s -1 . 

Figure 3 Spectra of two intergalactic regions in the cluster A2218, night 1. a) The inner, or 
“on,” aperture ( R/a = 0.94); b) the outer, or “off,” aperture (R/a = 5.33). Intensity 
is in units of 10 -18 ergcm -2 arcsec -2 A -1 s” 1 . 

Figure 4 Spectra of two intergalactic regions in the cluster A2256, night 1. a) The inner, or 
“on,” aperture ( R/a = 0.48); b) the outer, or “off,” aperture ( R/a, = 2.96). Intensity 
is in units of 10 -18 erg cm -2 arcsec -2 A -1 s -1 . 

Figure 5 Spectra of two intergalactic regions in the cluster A1413, night 2. a) The inner, or 
“on,” aperture ( R/a = 0.65); b) the outer, or “off,” aperture ( R/a = 2.94). Intensity 
is in units of 10 -18 erg cm -2 arcsec -2 A -1 s -1 . 

Figure 6 Spectra of two intergalactic regions in the cluster A2218, night 2. a) The inner, or 
“on,” aperture ( R/a = 0.94); b) the outer, or “off,” aperture ( R/a = 5.33). Intensity 
is in units of 10 -18 erg cm -2 arcsec -2 A -1 s -1 . 

Figure 7 The on-off spectra corresponding to Figs. 2-6. In each of these the spectrum in panel 
b has been subtracted from the spectrum in panel a. Cosmic-ray hits and poorly 
subtracted night sky lines have not been removed, a) A1413, night 1, b) A2218, night 
1, c) A2256, night 1, d) A1413, night 2, e) A2218, night 2. Intensity is in units of 
10 -18 erg cm -2 arcsec -2 A -1 s -1 . 

Figure 8 Three on-off spectra with an axion decay line obeying eq. (14) (with C = 1) artificially 
introduced into each spectrum. Intensity is in units of 10 * 8 erg cm arcsec A s 

a) A2256, night 1 with a 3.2 eV axion decay line added (indicated by arrow). This 
spectrum constitutes the low mass limit in this axion search (compare with Fig. (c). 

b) A2218, night 1 with a 3.5 eV line added (indicated by arrow) (compare with Fig. 
7b). c) A2218, night 2 with a 6.5 eV line added (compare with Fig. 7e and Fig. 6). 

Figure 9 The on-off spectrum of A2256, night 1 with obvious cosmic-ray hits and poorly sub- 
tracted night sky lines removed (compare with Fig. 7c). 

Figure 10 The cross-correlation function, £(/), for two sets of clusters. The lag is in pixels and 
the correlation = h < 1. a) The night 1 data of A2218 and A2256 cross-correlated. 
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A correlation peak arising from intracluster line emission would be centered at a lag 
of 138 with a width of ~ 30 pixels (compare with Fig. 11a). b) The night 2 data of 
A1413 and A2218 cross-correlated. A correlation peak from intracluster line emission 
would reside at a lag of 27 pixels. The narrow peak at lag = 0 is due to residual night 
sky lines. 

Figure 11 £(/) for two sets of clusters, each set with an axion line artificially introduced into the 
spectra, a) A2218 and A2256, each with a 3.5 eV, £ = 1 axion line in it (see Fig. 8b). 
The lag is in pixels and correlation = h. The correlation peak is centered at lag = 
138 pixels and is highly significant (h/2cr^ = 11.9) (compare with Fig. 10a). b) £(/) 
for the night 2 data of A2218 and A1413 with a 6.0 eV, £ = 0.07 axion line added to 
each (see Fig. 8b). The peak is centered at lag = 27 pixels and is highly significant, 
V(> h) ~ 1.8 x 10 -5 ( h/2a ^ — 3.04) (compare with Fig. 10b). 

Figure 12 2-cr limits to any line flux arising in the five observed clusters of galaxies; see text for 
details, a) A1413, night 1, b) A2218, night 1, c) A2256, night 1, d) A1413, night 2, 
e) A2218 , night 2. The large region with negative flux in the night 2 data are due to 

problems with the dome flat in the blue, see text for details. Intensity is in units of 

10 -18 erg cm -2 arcsec -2 A -1 s -1 . 

Figure 13 The values of £ required to hide an axion line in the observed clusters of galaxies. The 
dashed line is from cross-correlating spectra with noiseless templates (see table 3). The 
solid line arises from the flux limits of Figs. 13a-e, see text for details. 

Figure 14 a) A night 1 spectrum from A2256 minus that of A1413. b) A night 1 spectrum of 

A2218 minus that of A2256. Intensity is in units of 10~ 18 erg cm -2 arcsec -2 A -1 s -1 . 
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Figure 2b 
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